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We investigate the problem of what evolutions an open quantum system described by a time-local 
Master equation can undergo with universal coherent controls. A series of conditions are given which 
exclude channels from being reachable by any unitary controls, assuming that the coupling to the 
environment is not being modified. These conditions primarily arise by dehning decay rates for the 
generator of the dynamics of the open system, and then showing that controlling the system can only 
make these rates more isotropic. This forms a series of constraints on the shape and non-unitality 
of allowed evolutions, as well as an expression for the time required to reach a given goal. We give 
numerical examples of the usefulness of these criteria, and explore some similarities they have with 
quantum thermodynamics. 


I. INTRODUCTION 

The ability to coherently control quantum dynamics 
has received considerable interest in the last few decades, 
both for its potential application in technology and 
the insight it provides to fundamental science [5]-[7]. It is 
therefore surprising that the question of what dynamics 
can be reached with unitary controls is poorly understood 
in open systems where interactions with the environment 
cannot be neglected HIE]. Current tools rely primarily 
on finding explicit solutions to control problems [TO] but, 
as these methods are typically computational expensive 
and do not always give definitive answers, it is often hard 
to decide if the failure to find a good solution is due to 
its non-existence or simply an insufficient search. Having 
clear, efficiently accessible criteria which rules out cer¬ 
tain evolutions avoids these problems and would help in 
the quest for improving the design and optimisation of 
devices for quantum computation, communication and 
sensing. 

The question of what dynamics can be reached with 
coherent controls in the case of finite dimensional noise¬ 
less systems has been answered with the use of algebraic 
tools from the theory of Lie groups [num. Attempts 
to generalise these methods to open systems have met 
considerable mathematical difficulties. The two principle 
results are an accessibility criterion [HHni, which de¬ 
scribes which directions can be explored for short times; 
and Lie wedges m [T7|, which provide a partial char¬ 
acterisation of the geometry of the reachable set but in 
general cannot be calculated exactly. Other approaches 
focus on finding approximate numerical solutions HQlllH], 
or explore the related question of state-controllability, 
where the interest is in the ability to map one state to 
another [TOHSU • another approach is to treat the sys¬ 
tem and environment on an equal footing and approach 
the infinite-dimensional problem directly [25j |26] . 

In this paper we investigate operator controllability 
from a geometric approach and characterise broad ranges 
of evolutions that a dissipative system cannot reach with 



FIG. 1. We illustrate the principle ideas of this paper by show¬ 
ing a cross section of state space at two different times, with 
the small arrows indicating the direction of flow induced by a 
Lindbladian (the generator of memoryless dissipative dynam¬ 
ics). The large arrows correspond to the available controls. 
We see that different parts of the space are contracting at dif¬ 
ferent rates; by rotating the system in time with Hamiltonian 
controls some of these decay rates can be averaged together. 


any unitary controls. In the case that there are ex¬ 
perimental constraints that limit which of these con¬ 
trols could be physically realised, there are additional 
constraints imposed. The validity of the conditions 
detailed here is not affected however, as reducing the 
set of allowed operations can not increase what can be 
achieved. Furthermore our approach does not need a de¬ 
tailed knowledge of the behaviour of the environment or 
of the controls applied, but does require that the action 
of the environment on the system remains unchanged. 
Ideas similar to those we are presenting were introduced 
in EH in the case of time-independent and unital quan¬ 
tum systems with a special focus on single qubits. 

In order to get an intuitive understanding for the prin¬ 
ciple idea behind our work, it is useful to consider dis¬ 
sipative quantum processes as a flow in state space, as 
illustrated in Fig The noise can act in a variety of 
ways on the state space, including rotating and shrink- 
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ing the space (corresponding to decay) in a potentially 
anisotropic way. Hamiltonian controls allow us to im¬ 
pose additional rotations on the system such that differ¬ 
ent parts of the state space feel different contraction rates 
at different times. This results in the ability to mix the 
decay rates together and leads to the overall evolution 
obeying some averaged rates. These cause the final state 
space, which represents the total evolution, to be more 
isotropic than in the absence of controls. 

Our main results relate directly to this, and state that 
a quantum operation cannot be reached if it has a more 
ordered structure than the noise acting on the system. 
After introducing the specific problem we are addressing 
and the relevant mathematics in section |n| we define the 
decay rates of an open system with this structure in mind 
in section IIT and show that the action of any coherent 
controls is to make these more uniform. We demonstrate 
that the sum of these rates is unaffected by control and 
thus obtain a strict condition for the times at which a 
target evolution can be reached. In addition, we obtain 
bounds on how the purifying power of noise can be en¬ 
hanced by control. The strength of these criteria is tested 
numerically in section [IV| for common examples of noise, 
and we show that, at least for small systems, the neces¬ 
sary conditions are strong and tight enough to provide a 
major restriction on what evolutions are possible in real¬ 
istic situations. Both the language and the mathematics 
used to describe these relations are reminiscent of ther¬ 
modynamics, a link which we explore in section [V| We 
conclude in |VI| with a summary of the results, a compar¬ 
ison of the methods of this paper with prior results on 
Lie wedges, and possible directions for future work. 


II. PROBLEM AND MATHEMATICAL 
BACKGROUND 

The aim of this paper is to obtain some general rules for 
which operations cannot be performed on coherently con¬ 
trolled dissipative quantum systems by formalising the 
intuition described in Figl^and applying it to a more gen¬ 
eral setting. To do this we first state the control problem 
formally. We consider Hamiltonian controls on a finite di¬ 
mensional system interacting with an environment, with 
the requirement that the reduced system obeys a time- 
local master equation 

j^p = G,ip) = G°,{p)-i[H„p] 

pT = Mt{Po) = (1) 

where p is a quantum state, Gt is the linear generator 
for the motion, T is the time-ordering operator, and Mt 
is the resulting dynamical map, the set of which (vary¬ 
ing over total times and controls) we aim to characterise. 
The generator is divided into an uncontrollable drift 
and a controllable Hamiltonian term Ht. The latter is a 
time-dependent control Hamiltonian chosen so as to gen¬ 
erate the desired dynamics, and we impose no restrictions 


on it beyond being Hermitian. G^ represents the intrinsic 
part of the dynamics, such as an internal energy splitting 
or an interaction with the environment. If we restrict it to 
be a Lindblad operator, which we will denote by Lt, then 
the allowed solutions are Markovian, completely-positive 
trace-preserving maps. Although the Lindbladian case is 
the most commonly used and the one with the clearest 
physical interpretation, the key results of this paper do 
not rely on the specific form of the Lindblad operator 
and hold for a more general generator which gives rise to 
non-Markovian dynamics. In such cases there can be sub¬ 
stantial additional restrictions beyond the ones presented 
here, since the dynamics induced by control Hamilto¬ 
nian and non-Markovian drift does not necessarily induce 
completely-positive dynamics, even if the uncontrolled 
evolution is completely-positive |28]. The validity of the 
conditions presented in this paper are, however, not im¬ 
paired by the additional intricacies of non-Markovian dy¬ 
namics provided that the dynamics are still describable 
by a time-local linear generator. 

An implicit assumption in Eq.Q is that the con¬ 
trol Hamiltonian does not affect the dissipative com¬ 
ponent of the generator. While there are cases where 
this approximation holds exactly [29], this is not always 
the case, and it is possible for controls to modify the 
decoherence induced by the environment [30l [3lj. In 
such circumstances a variety of different methods have 
been developed [32H36] . but these typically require de¬ 
tailed knowledge about the environment or additional 
assumptions about finite-dimensionality or bounded in¬ 
teractions. They may also call for experimentally dif¬ 
ficult regimes such as strong and/or rapidly oscillating 
control fields. For these reasons, it is highly desirable to 
explore what can be achieved with controls if the dissi¬ 
pative component of the dynamics is not modified, which 
is the regime studied in this paper. 

In order to make the picture introduced in Fig rig¬ 
orous, and to derive our results, it is necessary to in¬ 
troduce some mathematical concepts and notation. The 
starting point is to work in the generalised Bloch repre¬ 
sentation [37], where a quantum state p is represented 
as the real vector \p) = (xq, xi, X 2 ,..., Xc^ 2 _i)^, where 
Xi = Tr[cr^p] are the expectation values over an orthonor¬ 
mal set of traceless Hermitian matrices for i = 1,..., —1, 

(Jo = d is the dimension of the underlying 

Hilbert space. In this representation super-operators 
acting on states become matrices. We will denote dy¬ 
namical maps and generators in this representation by 
M and G respectively to distinguish them from their 
super-operator form. The spectral properties (eigen¬ 
values, singular values, trace and determinant) of the 
super-operators are given by those of their matrix rep¬ 
resentation. The dual of a super-operator M”*", defined 
according to Tr [/iM(p)] = Tr[M’*'(/i)p], has as its ma¬ 
trix representation the Hermitian conjugate of the matrix 
representation of the original super-operator, such that 
(Mt) = (M)t. 

Writing out the explicit form of M highlights some 
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of its properties. When the dynamical map is trace¬ 
preserving M is of the form 

f ^ 

vi Mil 

V2 M 2 I 

V. 

where the top row is fixed, 

M is Hermiticity preserving 
jority of physically sensible 
to the reduced matrix. This form has the advantage of 
explicitly separating the unital and non-unital part of 
the dynamics. Unitality refers to leaving the maximally 
mixed state unchanged; as this is the only state left in¬ 
variant by all Hamiltonians and it is the centre of rota¬ 
tions, this is an important property for control. The left 
hand column consisting of the elements Vi fully describes 
the non-unital part of the map and quantifies how much 
the maximally mixed state is translated by the dynam¬ 
ical map. For unital operations, such as unitary evo¬ 
lution, these vanish and the dynamical map reduces to 
M = li0M- The reduced matrix M thus describes solely 
the unital part of the evolution. These are partially de¬ 
coupled from the non-unital dynamics in the sense that 

MbMA = M^MA. meaning that the total unital dy¬ 
namics of a concatenation is given by the concatenation 
of the unital part of the individual super-operators. This 
can easily be seen by noting that the concatenation of 
super-operators is given by the product of their matrix 
representations. 

As unitality is a key property, it is also useful to have a 
measure of how unital or non-unital a dynamical map is. 
A convenient one is Tr ^ the purity of the state 

obtained by applying the map to the maximally mixed 
state. This value is maximised at one if the maximally 
mixed state is mapped to a pure state and is minimised 
to ^ if the map is unital. In a similar fashion the non- 
unitality can also be quantified by higher moments 
Tr , and the first d moments are linearly inde¬ 

pendent. These should all be seen as describing roughly 
the same physical quantity: the ability of the map to 
purify states. 

The matrix form of a trace-preserving generator G is 
identical to Eq. except that the entire top row van¬ 
ishes. Its unital and non-unital parts can be separated 
in a similar way and the same results on concatenation 
holds. Because of this, we have that = e— and, hence, 
that if M is the dynamical map generated by Gt then M 
is the one given by Gt- This means that contained inside 
every non-unital problem is a unital one with a dimen¬ 
sion 1 smaller, and any solution to the control problem 
of generating M must also solve M. It is possible for this 
reduced problem to be non-Markovian even if the origi¬ 
nal one is Markovian, but this does not affect the validity 
of the approach. 

As they are particularly important, we note the form 
that the super-operators of closed dynamics take in this 


M 

M, 


...\ 


/ 


( 2 ) 


all the elements are real if 
(which it is for the vast ma- 
cases), and the tilde refers 


representation. Unitary propagators, M(') = U •I/’*', be¬ 
come matrices in the defining representation of the ro¬ 
tation group li 0 SO{d^ — 1). Hamiltonian generators, 
G(') = —i[H’, •], are in the corresponding Lie algebra, 
0i 05o((i^ — 1), which consist of real antisymmetric (and 
therefore traceless) matrices. In the case of d = 2, uni¬ 
tary propagators form all such matrices (and like-wise for 
Hamiltonians), but in higher dimension they only form 
a subgroup. This means that in d = 2 only, the set of 
every unitary on the system corresponds to all possible 
rotations of the state vector \p). In higher dimensions 
however, there exist rotations of this vector which can¬ 
not be induced by any Hamiltonian controls on the sys¬ 
tem; this is equivalent to saying that not every vector in 
the generalised Bloch space is a valid quantum state [37] . 
This property is one of the fundamental reasons why the 
results of this paper are necessary conditions rather than 
a complete characterisation of the allowed dynamics; we 
can only rule out targets which cannot be reached by any 
rotations, whether these are physical controls or not. 

The intuitive picture described in the introduction re¬ 
lies on a notion of averaging a set (the decay rates) to 
obtain another. A natural way to describe this process is 
the majorization relation [38| which tests if one real set 
is more uniformly distributed than an other. A set a of 
real numbers is majorized by another such set 6 , written 
as a ^ 6, if and only if 


ai < hi 

ai + a\<hi + b\ (3) 

where | signifies that the elements of the set are sorted 
in decreasing order. Another way of stating this is that 
a is majorized by b if and only if an ordering of a can be 
obtained by a convex sum of different orderings of 6 . It 
is this property which makes it suitable to describe an 
averaging procedure. Another useful property is that it 
is conserved under scaling such that a < b also implies 
Xa < xb ioi all real (including negative) x. A point to 
note is that, unlike standard inequalities on the reals, 
majorization provides only a partial order. 


III. THEORETICAL RESULTS 

With this formalism we are now in a position to re¬ 
fine the intuition developed in Fig The picture was 
of a Lindbladian (or a more general drift) acting on the 
state space such that it flowed from one shape to an¬ 
other. By coherently controlling the system, the space 
can be rotated so that some of the decay rates are aver¬ 
aged together. For a Markovian two-level system where 
the Lindbladian has no Hamiltonian component, these 
decay rates are the eigenvalues of the Lindbladian, which 
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are always non-positive. In higher dimensions however 
these eigenvalues can be complex which gives rise to two 
problems as they may not faithfully quantify the contrac¬ 
tion of the space and it is not clear what averaging them 
would signify. The situation is exacerbated in both the 
non-unital and non-Markovian case where there are even 
fewer constraints on the spectrum of the generator. This 
points to the need for a different way of quantifying the 
decay rates of an open system than naively taking the 
spectrum of the drift. 

Worki ng thro ugh the mathematics in detail (as we do 
in section III A) shows that the correct decay rates to con¬ 
sider are the eigenvalues of the sum of the drift and its 
dual, or equivalently, of the Hermitian part of G. Their 
physical relevance is supported by two important prop¬ 
erties. Firstly, as the generators of rotations in G are 
anti-Hermitian (whether it corresponds to a Hamiltonian 
degree of freedom or not), these rates capture only the 
decay/growth component of the flow. Secondly they are 
real, so averaging them corresponds to a more uniform 
flow in a way which can be naturally defined using ma- 
jorization. These suggest that the eigenvalues of the Her- 
mitian part of the drift capture some of the key aspects 
of the controllability of the system, and their prominence 
in the criteria detailed below show that this is indeed 
the case. In a similar fashion it can be seen that the 
anisotropy of the dynamical map is described by its sin¬ 
gular values (as they are rotationally invariant and non¬ 
negative), which loosely correspond to the characteristic 
lengths of the final state space. 

This, together with majorization as described above, 
allows us to write the conditions which must all be satis¬ 
fied for a dynamical map M to be reachable by a system 
with drift G^. These are: an expression for the state 
space volume reached at a given time (independent of 
the controls). 


det(M)=e/o"Tr(G?)dt^ 


( 4 ) 


which gives rise to independent conditions forn = 2,..., d. 
In order to highlight the systematic similarities of the 
conditions, and stress that it is the Hermit ian part of the 
drift that matters (which in the Bloch representation re¬ 
moves any Hamiltonian or other rotational component), 

G^ can be replaced by in Eqs.(|4) and Q. 

Instead of seeing if the system can reach a target map, 
we can instead ask if a time-independent generator G^ 
and controls can approximately simulate another gener¬ 
ator G' arbitrarily well. In this case we obtain the con¬ 
dition 


A (G' + G'^) ^ A (GV G^^) (8) 

whose unital version 

A (o' + G'1') ^ A (g° + G°1') (9) 

also holds provided both generators are trace-preserving. 

We proceed to give a proof of these conditions, followed 
by a detailed discussion. 


A. Proofs 

Evolution time — To derive Eq. we begin by not¬ 
ing that the formal solution for M is given in terms of 
time-ordered matrix exponential, which can be expressed 
according to the Magnus expansion m 

M = 

— ^So G(ti)(iti+2 Jq dti dt2[G{ti),G^{t2)]-\-■ ■ ■ 

where all higher order terms in the series consist of nested 
commutators. As the determinant of a matrix exponen¬ 
tial is the exponential of the trace we can rewrite this 
as 


which extends prior work on the determinant of quan¬ 
tum channels [39]; a constraint on the anisotropy of the 
dynamical map. 



where A and cr refer to the set of eigenvalues and singular 
values of the super-operator respectively, and the log acts 
element-wise on the set; a unital version of this condition 
(for trace-preserving drifts). 



which was shown in m for a time-independent Lindbla- 
dian; and bounds on the maximal non-unitality that can 
be reached. 


Tr[M(il)”] < sup {Tr[p”] \ 3t: Tr[p^-^G‘i{p)] = O} 

P 

( 7 ) 


det (M) = G(ti)dti-\-^ So dti dt2[G(ti),G(t2)] + ...'^ 

= eTr[S^^ Git)dt] 

where we have used the fact that commutators are trace¬ 
less. As discussed previously, the control Hamiltonian 
appear in the equation of motion 0 as commutators, 
therefore their Bloch representation are also traceless giv¬ 
ing Tr[G(t)] = Tr[G^(t)]. As the trace and determinant 
of M and G are identical to those of M and G, this gives 
the desired expression 

det(M) = (12) 

In the case that the Magnus expansion does not converge 
(which may happen if \ \G{t)\\ 2 dt > it |40|), then the 
proof can be extended by splitting the propagator into 
sufficiently many terms 

M = Te'^Gi G{t)dtq-^J^^ G(t)dt 


(13) 
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such that the Magnus expansion converges for each term. 
Applying the same steps as before to each term and using 
the fact that the determinant of a product is the product 
of the determinants we arrive at 

det (M) = det (re-^‘n det 

= eftl Tr[G?]<it_ _g//= Tr[G?]dtg/o‘i Tr[G°]dt 
= Tr[G?]dt+...+/4^ Tr[G?]+/4i Tr[G?]dt} 

= efo Tr[G?]dt 


each term in the previous equation can be expressed for 
small 6t as 


logcr(e—= log 




log 


X yG°{t)+G°Ht))St+[G°{t),G°^t)]Sy 2 


= jA (iG°it) + G°\t))St + [G\t),G°\t)]5t^') 

(19) 


as before. 

Anisotropy of the dynamical map — The proofs of 
Eqs.([5][^ arise from two observations. Firstly, the evo¬ 
lution can always be decomposed into infinitesimal time- 
steps in a Trotter-like way, alternating between coherent 
and incoherent evolution. Secondly, the controls only af¬ 
fect the coherent steps which are all rotation matrices 
and so do not modify the singular values of the incoher¬ 
ent time steps, as singular values of a matrix depend only 
on the product of that matrix with its Hermitian adjoint. 
To prove Eq.([^ we expand the time-ordered exponential 
in terms of short time steps 

M = Te^o G(t)dt^ 

= lim (15) 

(5t^0 V J ^ ^ 

We now consider the singular values of both sides of the 
equation, denoted by the operator a. Specifically, we use 

m 


loga{AB) -< logcr(A) + loga{B) (16) 


for the majorization relation between the singular values 
of matrices and their products, where the log is under¬ 
stood as acting on each element in the set, and the sum 
on the right hand side acts on the elements of the sets 
ordered by magnitude. Generalising this to the case of 
multiple sums and applying it to Eq.(15), we obtain 


loga(M) lim I loga (e^°+ - (17) 

St^O I \ / 


As mentioned above, the coherent steps corresponds to 
rotation matrices and therefore do not affect the singular 
values. This allows the expression for the singular values 
to be simplified to 


loga(M) ^lim^ |log[t7(e—+ ... -1- log[t7(e—| . 

(18) 


We recall that singular values are obtained by a (A) = 
A ( V AAt), where A signifies the eigenvalues. From this. 


where higher order terms can be calculated using the 
Baker-Campbell-Hausdorff formula. The first term is of 
order 6t and, as the number of terms in Eq.(18) is 
it contributes to the integral in the limit St ^ 0 while 
all the higher order terms vanish. This gives as the final 
expression 


( 20 , 

which is independent of the representation used and 
so holds for the super-operators themselves. The 
proof for condition (|^, the unital version of this 
for trace-preserving generators, follows immediately 
from the fact that M = Texp{ G(t)dt} implies 

M = TgxpIJq^ G(t)(it}, as was discussed in section [ll} 
A different proof of this latter unital result was shown m 
[27] for time-independent Lindbladians only, and relied 
on similar mathematical ideas. 


Maximal non-unitality — To prove the non-unitality 
bounds, Eq. 0 . we begin with the formal expression and 
then derive an easily evaluable bound for it. To do this 
we are required to make the additional assumption that 
Gt is continuous. The maximal non-unitality of an open 
system is quantified by 

Tr[M(il)”] < supTr[p”(t,if,)] (21) 

t,H-r 

where p{t,Hr) = Te>^o 


and the supremum is over all possible evolution times 
and all possible controls. It is sufficient to consider only 
Hr which are defined for r G [0,oc). Eq.(21) appears 
as difficult to calculate as solving the control problem, 
and is therefore of limited use. However an upper bound 
for it can be found more readily. To do this we note a 
property that the supremum must satisfy as a function 
of t for any Hr. From this, we reformulate the constraint 
that p has been evolved from the maximally mixed state 
into one which is easier to work with. 
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For a given continuous finding the supremum of 
Eq.(21) reduces to finding the supremum of a scalar func¬ 


tion which is bounded between ^ and 1 and differen- 

d 

tiable everywhere. There are several cases in which this 
could happen. Firstly, the supremum being reached at 
t = 0 can be immediately excluded as the function is a 
minimum at that point. Secondly, the supremum being 
reached for some finite time leads to J^Tr[p’^(t)] = 0 at 
that point in time. Lastly, if the supremum is not reached 
for finite t then either: it is reached in the limit t ^ oc 
and so the derivative also goes to 0 in this limit (due to 
the function being bounded from above), or the limit is 
undefined because the function does not converge. In the 
latter case there are many local maxima which form a se¬ 
ries, the supremum of which gives the supremum of the 
original function. As the gradient of each these maxima 
is 0, the largest value that can be reached by the function 
also occurs when the gradient vanishes. Hence, we have 
that a necessary condition for the supremum of Eq.(21) 


iTr[p^]=0, 

Tr[p"-i G°ip)] - iTr [p^-^ [H^ p]] = 0, 
Tr[p"-iGO(p)]-iTr[[ift,p"]]=0, 
Tr[p"-iG0(p)] =0. (22) 


If this condition is satisfied for infinitesimal 5t then, by 
concatenation, it holds for all time T and can effec¬ 
tively simulate G'. We use the term ‘effectively’ to em¬ 
phasise that although the generator G' cannot be reached 
exactly, it is possibly to follow a trajectory in state space 
which is arbitrarily close to the one generated by it. Un¬ 
der these conditions Eq.(24) simplifies to 


loga(e^ '5*) |A (G° + 6t, 


logA(e®'^*e^''^*)5 ^ 




logA 


(25) 


A {C + G'^) ^ A (G° + G°^) . 


The unital version of this relation also holds provided 
both generators are trace-preserving for the reasons dis¬ 
cussed in section im We note that this condition implies 
that Tr[G'] = Tr[^] is also required. 

There are cases where the time it takes to simulate 
the dynamics is not of concern, which corresponds to the 
traces of G and G' not being equal. In such cases, the 
condition above can be relaxed to 


Since the gradient as calculated above depends solely on 
G^ and not on the controls, we can relax the condition on 
Gt being continuous to G^ being continuous. Instead of 
calculating the supremum over all controls, we can com¬ 
pute it over all states which satisfy this condition. This 
allows us to place a bound on the maximum non-unitality 
reachable by a dynamic system which, as desired, does 
not require any propagators to be calculated 

Tr[M(il)^] < sup {Tr[p^] \ 3t: Tr[p^-^G?(p)] = O} . 

(23) 

This provides up to d constraints (including the trivial 
case for n = 1), as higher moments of p are not indepen¬ 
dent. It is interesting to note that similar expressions 
were also arrived at in a different control problem, that 
of finding the possible steady states of a driven open 
system [24] . 

Generator anisotropy — The criterion for generator 
anisotropy (applicable only to time-independent genera¬ 
tors) is similar to the one for the anisotropy of the dynam¬ 
ical map, with the important difference that the target is 
a flow in state space which we desire to achieve contin¬ 
uously in time, rather than a snapshot of the evolution 
at a single instance. The derivation for this condition 
begins with Eq. (§, where we replace the target M by 
e^ ^, and limit ourselves to the system drift being time- 
independent, such that our starting point is 

log[<,(e°'’-)| ^ IP vr. (24) 


2Trb']A ^ 2Tr[Gi’]'^ (^° + (^6) 

by a rescaling of time in Eq.([24|). This comes by 
replacing T with on the left hand side that that 

equation. This is a relaxation of Eq.([^, as it holds even 
if Tr[G'] 7 ^ Tr[G^]. 

Unital qubit Lindbladians — Eor the case of qubits 
undergoing unital Lindbladian dynamics, Eq.([^ can be 
simplified to m 

A (V) ^ A (L°) (27) 

as the dissipative part of the Bloch representation of 
any unital qubit Lindbladian is symmetric [41] and drift 
Hamiltonians are not of interest to us (they can be can¬ 
celled out by controls). This criterion is also sufficient. 
To prove it, we provide an explicit way to reach L' using 
a drift Lq and unrestricted Hamiltonian controls. Eor 
simplicity, we pick the time scale of the target dynamics 
such that Tr[L'] = Tr[Zo] = 1. By using the singular 
value decomposition, the target can be expressed as 

0 0 \ 

M = = UDV = U\ 0 0 V (28) 

\ 0 0 

where the —ni are the eigenvalues of L , and U and V 
are elements of 0(3). Eurthermore, as M has a positive 
determinant and the diagonal block is positive, we can 
pick U and V to have determinant +1, thereby restricting 
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them to SO(3). In a similar way, we can express the free 
evolution of the system for time t as 

0 0 \ 

eLot _ WF{t)W'^ = W \ 0 0 (29) 

\ 0 0 


where the —jii are the eigenvalues of Lq. Using the same 
argument as above, W can be chosen to be in SO(3). 
Controls on the system allow the implementation of any 
R = e— which, as we noted previously, corresponds to 
any matrix in SO{3). The control scheme to reach the 
target map corresponds to alternating free evolution and 
instantaneous controls as 


UDV = RiWF{ti)W'^R2...WF{tn)W'^Rn+i (30) 


We pick Ri = UWt, = WV and relabel RjW = 
R'j which can always be done due to the group structure. 
Next we pick the R' to be permutation matrices (which 
all he in SO(3)) such that Eq.(30) consists solely of diag¬ 
onal matrices where every term is an exponential. This 
lets us express the previous matrix equation in the simple 
form 


I ^2 1 = I 1 U + I /^3 1 ^2 + ••• + I /^2 1 ^6: (31) 

V^3/ V^3/ V^2/ V/il/ 

where, we recall from the way we picked the scale of Lq, 
that = 1. This control scheme thus allows us to 

reach any M = e- ^ where the eigenvalues of L are a 

convex combination of those of L ; which is equivalent to 
saying that they are majorized by them [38]. This means 
that A (I/') ^ A (l/^) is a sufficient, as well as necessary, 
condition for reachability in unital qubit systems with 
unconstrained Hamiltonian control. 


B. Discussion 

Evolution time — The first criterion, Eq. is 
an equality which appeared in [39] for the Markovian 
and time-independent case. Extending it to the time- 
dependent case gives it an important use in control the¬ 
ory: it states that a target map may only be reached by a 
dynamical system for the times which satisfy Eq.Q. To 
understand why this is the case we note that the modu¬ 
lus of the determinant of M is the volume occupied by 
its image and the trace of Gt is the rate at which this 
state space is growing (this will be non-positive, unless 
the system is non-Markovian, leading to a contraction of 
the space). The interpretation of this result is thus that 
the total rate at which volume is lost in state space is in¬ 
dependent of the Hamiltonian controls. More insight can 
be gained by noting that the trace of the drift is always 
positive in the physically sensible case of the generator 


being Hermiticity preserving. This means that the evolu¬ 
tion of such a system can only reach maps with positive 
determinant usms]. As this is not the case for every 
completely-positive trace-preserving map, this condition 
allows us to immediately rule out large sections of the 
space as unreachable for a broad class of dynamical sys¬ 
tems. 

If we further restrict the drift to be a Lindbladian 
(with a non-vanishing dissipative part) at all times, then 
the trace is always negative, signifying that a target 
map can only be reached at a single instant in time 
(if at all), and that this time can be easily calculated 
as it is independent of the controls. If the drift is 
non-Markovian (where the interplay between memory 
effects and controls has received much recent attention 
[42l [43|) the trace of the drift can be positive for certain 
times, leading to revivals in the determinant. Indeed, 
this has already been suggested as an indicator of 
non-Markovianity [44|. Eor our purposes, this particular 
feature leads to the possibility of there being several 
solutions to Eq. @ for a given target map. In both 
cases, the precise information about required evolution 
time given by this condition is in stark contrast to 
the case of closed systems, where in general very lit¬ 
tle is known about the time required to reach a target 
without explicitly solving for the evolution of the system. 

Anisotropy of the dynamieal map — The conditions 
of Eqs. § and § are a refinement of the intuition that 
controlling the system allows us to average the decay 
rates of the drift together. This is most easily seen by 
considering the unital case. The left hand side of the 
relation are the singular values of the dynamical map 
which, in the same way as the determinant is the volume 
in state space, are the characteristic lengths of the final 
state space. Thus, while Eq. 0 determines the volume 
reached, Eq. 0 provides a constraint on the anisotropy 
of the dynamical maps that can be reached. 

The non-unital majorization relation has broadly the 
same interpretation, although the overall shift caused 
by the non-unitality manifests itself in the decay rates 
and singular values in a complex way. Indeed, one of the 
eigenvalues of + G^^ will typically be positive in the 
non-unital case, even if the drift is Markovian. While the 
idea of a positive decay rate in a Markovian system may 
appear counterintuitive, it only signifies that some states 
become purer under such a Lindbladian. This is most 
easily seen by considering the Bloch sphere: negative 
eigenvalues correspond to states moving towards the 
maximally mixed state, but if the system is decaying 
to the state |0) then there is also a dynamic evolution 
away from the centre towards a pure state on the 
surface of the sphere. Thus, despite their similar form, 
Eqs.([5]) and 0 give very different results and there 
are many dynamical maps that satisfy one but not 
the other for a given drift (as we will show in Eig [^. 
It is also worth noting that Eq. 0 is recovered, up to 
a modulus, by the last term in the majorization relations. 



Maximal non-unitality — This last condition on dy¬ 
namical maps, Eq.Q, is conceptually very different from 
the others. Rather than restricting the shape of the dy¬ 
namical map, it provides a series of constraints on how 
much the maximally mixed state can be displaced, cor¬ 
responding to where the centre of the image of the dy¬ 
namical map lies in the state space. Although the right 
hand side of Eq. 0 is independent of controls, the max¬ 
imisation over all states (and over all t if the drift is 
time-dependent) makes this criterion somewhat harder 
to evaluate in higher dimensions. It is worth noting that 
the constraint is less strict than G? (p) = O 5 which means 
that the non-unitality is not bounded by the fixed points 
of the drift. The interpretation of this criterion is there¬ 
fore that it is possible to increase the ability of noise to 
purify states by using controls, but only up to the limits 
given. 

A physical example of this is a three-level system in a 
A configuration (such as in Eig|^, with the excited state 
decaying into the two low level states. In the absence 
of controls, the system has some non-unitality as the 
maximally mixed state over the three levels will decay 
to a mixed state over only two levels. With the use 
of controls, however, the population can be coherently 
transferred back from one of the two ground states to 
the excited state where it will once again decay. Doing 
this many times results in the total population being 
transferred to the other ground state and the total 
action of the dynamics is to map everything to a single 
pure state. Thus, this dynamical map induced by a 
specific set of controls has maximal non-unitality. This 
is the principle behind optical pumping and shows that 
non-unitality can be increased with controls. If the two 
lower levels had some decay between them, however, 
this scheme may not work perfectly and Eq. 0 provides 
bounds for how well it can be done. 

Generator anisotropy — Instead of investigating if 
the system can reach a target map, in Eq.(|^ we con¬ 
sider if it can be made to approximate a different drift 
continuously in time. To do this we limit ourselves to 
time-independent and see if it can give rise to evolu¬ 
tions arbitrarily close to M = e^ ^ for all t. If it can, we 
say that can effectively simulate G' as it can replicate 
the same dynamics arbitrarily well in a time continuous 
fashion. Necessary but not sufficient criteria to do this 
are given by Eqs.(|^ and 0. These are stricter than the 
anisotropy conditions on dynamical maps; it imposes not 
only a target map but the whole trajectory in time to it. 
That this can be done at all is at first hand surprising, as 
the only generators which can be reached exactly from a 
given drift are given precisely by the drift plus all possi¬ 
ble controls. However, by quickly rotating the system it 
is possible to get arbitrarily close to the required trajec¬ 
tory by winding tightly around the desired path without 
ever moving exactly along it. 

In the case of the generator G^ being a unital 


qubit Lindbladian, Eq. 0 can be simplified further to 
A (I/') ^ ^ and is suffieient [27]. Eurthermore the 
eigenvalues of such unital qubit Lindbladians, A(I/), 
are constrained by complete positivity [45]. Erom this 
it is straightforward to show that \{L) = — (|,^,0) 
majorizes the spectrum of all other such Lindbladians. 
Hence, one with such a spectrum, such as dephasing 
!/(•) = —[(Jz^[(Jzi’]]i is universal and it can simu¬ 
late all other unital qubit Lindbladians. Conversely, 
the completely depolarising channel with eigenvalues 
A(I/) = — (|, |) is majorized by all other Lindbladians 
and is therefore at the bottom of the hierarchy defined 
by majorization. 


IV. NUMERICAL RESULTS 

The conditions detailed above are almost all necessary 
but not sufficient, so the question of how tight they are 
is important. We carried out numerical simulations to 
quantify this in two ways: firstly by how often the cri¬ 
teria forbid a target from being reached, and secondly 
by how often a target can be reached when it is not ex¬ 
cluded. The issue with doing this is that, for the very 
reason that the criteria derived in this paper are useful, 
it is computationally very difficult to test if a dynamical 
map can be reached with a given drift and controls. The 
only definite method requires simulation and optimisa¬ 
tion of the control problem. The size of the simulation 
itself scales as and the cost of optimising the control 
pulses scales far worse uni. Nevertheless, we obtained 
results for a class of common non-unital qubit and qutrit 
Lindbladians. 

Eor the qubit system the Lindbladian we consider is 
generalised amplitude damping, a ubiquitous type of 
noise, corresponding to a qubit which can exchange an ex¬ 
citation with a bath at finite temperature m- This can 
be thought of as a spin which has a finite rate for transi¬ 
tioning from the excited to the ground state and from the 
ground to the excited state, where the ratio between the 
two is a function of temperature. The temperature deter¬ 
mines the non-unitality of the noise: at zero temperature 
the steady state is the pure ground state, while at infinite 
temperature it is the maximally mixed state. As the non- 
unitality of the Lindbladian is an important aspect of its 
controllability, we use the purity of the steady state as a 
parametrisation of temperature. We analysed the ratio 
of randomly generated time-dependent Markovian maps 
[48] which could be reached numerically, and whether 
they satisfied the criteria Eqs.([^-Q, for different values 
of the non-unitality of the drift. Eq.Q was used in de¬ 
ciding the evolution time for which we attempted to find 
solutions of the control problem. 

The results for such a Lindbladian at different temper¬ 
atures are shown in Eig Taken together, the criteria 
state that over 90% of the space is unreachable at each 
temperature considered, showing that the conditions are 
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FIG. 2. The graphs characterise the strength of Eqs.(5|7) for 
generalised amplitude damping Lindbladians as a function of 
the purity of the drift steady state - itself a function of the 
bath temperature. The top graph shows what fraction of the 
randomly generated channels satished all of Eqs.( [5][7t sind, 
out of those, how many could be reached with a numerical 
optimisation package |46] to within a distinguishability (given 
by the diamond norm) of at least 0.1%. The bottom graph 
shows what fraction of the same channels are ruled out by 
each of the conditions individually. 


useful as they cut out the large majority of dynamical 
maps as impossible to achieve. The insufficiency of the 
criteria manifests itself in that - at some temperatures - 
only 10% of those not ruled out can be reached. This 
number, however, approaches 100 % in the unital case, 
which is expected as we know that the majorization con¬ 
dition is sufficient in unital qubit systems. That this 
figure rises again for highly non-unital, low temperature 
baths shows that the criteria are increasingly useful in 
this limit too. It is also interesting to note that the rel¬ 
ative importance of the different conditions varies with 
temperature: when the noise has a pure fixed point the 
unitality criterion provides no information and the unital 
majorization criterion is the most restrictive, while their 
importance is reversed when the fixed point is maximally 
mixed. 

The second example is another common type of noise, 
sketched in Fig is a qutrit in a non-symmetric A con¬ 
figuration where the top level decays to the two lower 


levels according to the Lindbladian 

Lip) =J2Li {LipL\ - HLlLi , P}) (32) 

i=l,2 

where Li = \i) (3| [49]. We focus on how the controls and 
the skew, 71 / 72 , can influence the asymmetry of the final 
evolution. To do this we picked a drift with a fixed skew, 
and investigated how close the system could get to maps 
generated by a similar drift but with a different skew. 
Fig shows that those which were as or more symmetric 
(a skew closer to 1 ) as the drift could be reached with 
a very high fidelity, and increasingly poorly those which 
were less symmetric. This is in excellent agreement with 
the majorization criteria as plotted. The tightness of 
the necessary conditions in this scenario show how useful 
they are in cases where there is a clear measure of non¬ 
uniformity, demonstrating that we can use controls to go 
from a highly ordered evolution to a less ordered one, but 
not the other way around. 

This result may at first hand appear to contradict the 
conclusion arrived at in section IIIB| when this example 
was also discussed on maximal non-unitality in the con¬ 
text of optical pumping. There we said that a qutrit in 
a A configuration could have a pure fixed point regard¬ 
less of the ratio of the decay rates; while here we stress 
that the skew cannot be increased. The resolution of this 
problem is that although the fixed point of the dynamics 
can be chosen independently of the skew, this only de¬ 
termines the evolution at t ^ 00 , at all other times the 
state space occupies a finite volume and the shape of this 
volume is what is constrained by the skew. 

We expect that in general the tightness of the criteria 
would decrease as the system size increases. The main 
argument for this is that the number of rotations in 
the state space that do not correspond to Hamiltonian 
degrees of freedom grows with dimension. That the 
conditions were tighter in the qutrit than the qubit 
example, however, suggests that there are specific cases 
where they remain an excellent approximation to the 
allowed operations. Even when less tight, the potential 
usefulness of the conditions may be greater for larger 
systems as it is substantially harder to simulate these 
and therefore to learn about their controllability via 
other means. 


|3) 


| 1 > 


71 ^ ^72 


| 2 > 


EIG. 3. Sketch of a qutrit in a A configuration, where the 
Lindbladian causes the excited level to decay to the two 
ground state at different rates. The skew, 71 / 72 , characterises 
the asymmetry in the system. 
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Skew (Y^lY2) 


FIG. 4. The graph shows how the dynamical map anisotropy 
conditions relates to the controllability of the non-symmetric 
A system outlined in figure The system drift had a skew 
of 10 and we attempted to reach maps generated by the same 
Lindbladian but with skews between 1 and 20. Plotted is 
the minimal distance (given by the diamond norm) to targets 
with different skews that could be reached with a numerical 
optimisation package [46]. The thick red line at 10.3 is the 
non-reachable boundary given by the majorization criteria - 
everything to the right of it is excluded - in excellent agree¬ 
ment with the numerical results. The slight bump just be¬ 
low 10 is due to the difficulty of numerically finding solutions 
which require rapidly oscillating control Hamiltonians. 


V. RELATIONS TO THERMODYNAMICS 

The focus of this paper has been the majorization con¬ 
ditions and which provide limits on the oper¬ 
ations that can be reached on an open quantum system. 
It is interesting to note that majorization also plays a key 
role on the related question of state, rather than opera¬ 
tor, controllability. A central result there is that a state 
a can be reached from a state p by a unital completely- 
positive trace-preserving map if and only if A(cr) ^ A(p) 
m This is a stricter form of the Second Law of Thermo¬ 
dynamics as the majorization relation A(cr) ^ A(p) im¬ 
poses the constraint S{a) > S{p) on the von-Neumman 
entropy m of the two states, but the converse is not al¬ 
ways true. This means that, under unital evolution, the 
eigenvalues of a state cannot become more ordered and 
therefore the entropy cannot decrease. 

The natural extension of this to the results on gen¬ 
erator controllability, Eqs. (§ and (|^, is that it lifts a 
form of the Second Law from applying to states to super¬ 
operators. This is a restatement of what we have shown: 
that generators of dissipative dynamics cannot become 
more ordered by the presence of coherent controls. These 
relations imply that the process is irreversible, controls 
can be used to make an existing generator of noise ar¬ 
bitrarily close to a different generator G', but the reverse 
cannot be done even approximately. This is a surprising 


result as the controls themselves are fully reversible as 
they are coherent. The rise of irreversibility from purely 
reversible pieces is a long standing puzzle of quantum me¬ 
chanics and a key aspect of the Second Law, the criteria 
we have developed here shows that it applies to super¬ 
operators as well as states. 

A more explicit link between the present results and 
thermodynamics can be found by considering the rate of 
change of the entropy of a system as it undergoes unital 
evolution. An expression for this in terms of the spectral 
properties of the channel is given in [50] which can be 
easily modified for Markovian channels, using Eq.([^, to 
give 

(33) 

where p{t) is a state evolving under a unital Lindbladian 
(possibly under the presence of controls), || • ||2 is the L 2 
norm, and Ai is the smallest (in magnitude) eigenvalue of 
the unital part of the Hermitian part of the Lindbladian. 
While the left hand side is any channel generated by a 
Lindbladian and control, the lower bound is the smallest 
decay rate of the Lindbladian and independent of the 
controls. This shows that the minimal rate of entropy 
production cannot be lowered. The physical picture is 
that rotating a system as it decays cannot increase how 
well states are shielded from the production of entropy. 


VI. CONCLUSION 

The principle idea behind this work is that the decay 
rates of the generator of an open system, the eigenval¬ 
ues of the sum of the generator and its dual, provides 
limitations as to the operations the system can achieve 
with coherent controls, resulting in Eqs. The decay 

rates can be made more isotropic by coherent controls, 
corresponding to the rates being averaged out by rota¬ 
tions, but the total rate of decay cannot be changed and 
it is not possible to create a more ordered structure or 
to increase the non-unitality beyond a given limit. These 
hold for a range of open quantum systems - going be¬ 
yond Markovian ones - within some assumptions which 
are discussed earlier. 

These assumptions are shared with existing work on 
the controllability of Lindbladians based on Lie wedges 
m [i7]. The Lie wedge provides a sufficient but not 
necessary condition for controllability (as the semigroup 
closure still needs to be taken), while the criteria of this 
paper are necessary but not sufficient. Taken together, 
they allow us to approximate the reachable set from both 
sides. Our approach gives results which are easier to use 
and can be calculated numerically, while in many cases 
there are no known methods to determine the exact Lie 
wedge, especially in the non-unital case. It also has the 
considerable advantage of allowing drifts which are non- 
Mar kovian and time-dependent. However, the method 
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used here does not enable us to see what effect reducing 
the allowed set of controls has. In the simplest case of 
unital qubit Lindbladians, we saw that majorization was 
sufficient if we had unrestricted Hamiltonian controls; a 
result which can also be obtained from Lie wedges, show¬ 
ing the consistency of the two methods. 

The results we have presented are a useful tool in the 
quest for designing quantum systems to achieve desired 
non-unit ary tasks, as they rule out some dynamics as 
impossible without the high cost of simulation and opti¬ 
misation. Two examples of the use of the criteria were 
investigated numerically highlighting that, although 
they are necessary but not sufficient, they still give 
a practical approximation to the allowed operations. 
Due to the partial order induced by majorization, it 
shows that some types of noise are “more useful” than 
others as they can be used to replicate all the same 
evolutions, in addition to others. Directions to develop 
this further include applying it to different types of 


systems, such as Gaussian ones. Another avenue would 
be to investigate the effect of reducing the allowed set of 
controls, for example by considering the system as multi¬ 
partite where only controls local to the sub-systems 
are possible. Further investigations into the link with 
thermodynamics may also prove fruitful. 
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